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Abstract 

We report here observations of the active galactic nucleus Mrk 501, at energies above 250 GeV carried out 
with the CAT atmospheric imaging telescope from March 1997 to Autumn 1998. This source was in a high 
Q"^ ■ state of activity at several different wavelengths in 1997, and the observed flux at TeV energies has been seen 

to change by a factor of ~ 20 from from 1995 and 1996 fluxes. CAT observations also indicate a curved 
C I spectrum at TeV energies, and a correlation between the gamma-ray intensity and the spectral hardness. The 

temporal variability and the TeV spectral properties are examined. 
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^ ■ 1 Introduction: 

The active galactic nucleus (AGN) Mrk 501 is one of the closest (z=0.034) BL Lacertae objects. It was 
discovered at TeV energies by the Whipple Observatory (Quinn et al., 1996), and confirmed by HEGRA 
■^ ; (Bradbury et al., 1997) and the CAT collaboration (Punch et al., 1997). 

O ' From March to October 1997, VHE observations of Mrk 501 revealed extreme variability, with a measured 

Q ■ flux up to 8 times the Crab Nebula flux. Several independent Cerenkov telescopes have confirmed this dramatic 

Q\ ', activity (Protheroe et al., 1998). Based on observations carried out in 1997, we present a detailed study of 

^ I spectral properties. Thanks to the large variability of the source during this period, the spectral hardness was 

studied as a function of the intensity of the source and spectra were derived for different intensity levels. 
Mrk 501 has been regularly observed by the CAT telescope since March 1997 and the resulting light curve 
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O I is presented. 
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^ ; 2 Cat Cerenkov Imaging Telescope and data analysis 

Located on the Themis solar plant in southern France (2°E, 42°N, 1650 m a.s.L), the CAT (Cerenkov Array 
at Themis) imaging telescope has been described in detail in Barrau et al. (1998). It achieves a low threshold 
(~ 250 GeV) despite its small reflector area (17.7 m^) by taking full advantage of the rapidity of the Cerenkov 
C^ [ pulse with a near-isochronous mirror, fast phototubes, fast trigger and readout electronics. Moreover, its very 

high definition camera (546 pixels with an angular size of 0. 12°) allows an accurate analysis of the longitudinal 
and lateral light profile of the shower image, as discussed in Le Bohec et al. (1998), giving a good separation 
of gamma-ray showers from hadronic ones by means of a x^-like variable and of the pointing angle a. 

Cuts used in the standard CAT analysis are the following: P (x^) > 0.35, a < 6°, Qtot (the total charge 
in the image) > 30 7e ; an additional cut requiring at least 3 je in the fourth-brightest pixel ensures uniform 
trigger conditions. These cuts yield a background rejection factor of about 200, with a 7-ray efficiency better 
than 40%. 

Observations are done in an ON/OFF mode: the source is tracked (ON) for ~ 30 minutes and, to estimate 
the background, the telescope tracks a control region offset in right ascension on the source. Typically, one 
OFF-region is observed for every two ON-regions. Data are further selected for good sky-conditions and stable 
detector operation. 

3 Source Variability and Light Curve 

Here, we include data from Spring 1997 to Autumn 1998, representing a total of ~ 100 hours of observa- 
tion. The data-set has been further restricted to zenith angles < 45°, for which the detector calibration studies 
have been completed. Figure 1 shows the Mrk 501 nightly integral flux levels above 250 GeV. 

In 1997, VHE 7-ray emission was very strong and significant variability was observed in the nightly 



averages. The highest flare, recorded on 
April le*"^, 1997 reached about 8 times the 
Crab level. On the other hand in 1998, the 
integral flux appears to have been fairly 
constant and low (~ ^ of Crab nebula 
rate). 

While variability at a nightly scale is di- 
rectly seen in Figure 1, the search for 
intra-night variability has been studied in 
detail in Renault C, Renaud N. and Henri 
G. (1998). No significant short-term vari- 
ability was found. 
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Figure 1: Mrk 501 nightly integral flux levels above 250 GeV, expressed in 
units of 10~^^ cm~'^s~^, for observations in 1997. The average spectral shape 
as derived in Sect. 4 has been assumed. 



4 Differential Energy Spectra 

The effective detection area and the energy resolution function are derived from simulations as functions 
of energy and zenith angle. These simulations take account of the detector response and are calibrated on the 
basis of the large statistics of gamma-ray excess events from runs on Mrk 501 with a high signal-to-noise ratio, 
and of Cerenkov rings induced by muons. 

The differential energy spectra are obtained, for a given hypothesis on the spectral shape, by a 

E(TeV) 
1 10 



maximum-likelihood fit taking into 
account the effective area and the 
energy-resolution function of the tele- 
scope. Two possible spectral shapes 
have been investigated: a pure power 
law (Wo) and a curved shape (Hi) 
defined by ^q^tcV^'''"''"'''^"^^' ^ug- 
gested by general considerations on 
7-ray emission from blazars. The 
likelihood ratio of the two hypothe- 
ses, defined as A = — 2 x log ^W^. 

gives an estimate of the relevance of 
7^0 with respect to TYi. 
For this spectral analysis, the data- 
set was taken from March to Octo- 
ber 1997. During this period, three 
different activity states of Mrk 501 
have been defined: HF for high- 
intensity flares with an integral flux 
(> 50 X 10^1^ cm-^s-i), LF for 
low-intensity runs with an integral 
flux (< 12 X IQ-i^ cm^^g-i) and 
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Figure 2: vF„ or E'^ d<I>/d_B spectra for the HF, MF, and LF states. 
Individual intensities per bin are only indicative. 



MF for mid-intensity flares with an integral flux between HF and LF. The average over the complete data-set 
(AV) has also been considered. Results for each data-subset are given in Table |I|. Only statistical errors are 
quoted in this table. 

The resulting differential energy spectra of Mrk 501 are shown for energies from 330 GeV up to 10 TeV in 
Figure 2. The likelihood ratio, A, behaves like a x^ with one degree of freedom. Thus, from the values of A, 
the probabilities of falsely accepting TLi, over TIq, are 7 x 10~^ and 7 x 10^^^ respectively, if the power-law 
hypothesis were true. The weak intensity of the source in the LF state may explain the lower significance of 
the curvature term [3 in this state. The VHE peak emission of Mrk 501 clearly takes place in the range of 
several hundred GeV. The corresponding peak energies, -Emax. in the spectral energy distribution (£'^d(/)/dE) 
have been derived from the fitted values of a and j3 (Table 1). 

There is some indication of an increase in E'max with source intensity, which is equivalently shown by 
the variation of the fitted values of q, corresponding to a possible correlation between intensity and E'max- 
The intensity-hardness correlations, discussed in the next section, support this effect. Moreover, a correlated 
variability in X-ray and 7-ray have been observed (Djannati-Atai et al., 1999). So, together with the power 
deficit in the GeV range (Samuelson et al., 1998), this strongly suggests a two-component emission spectrum 
for Mrk 501 (Fossati et al., 1998). 

Table 1: The differential flux is given in units of 10^^^ cm^^s^^TeV^^. In the hypothesis of a curved 

shape (T^i), it is parametrized as d(^/d£'TeV = ^o^TeV °^^° ^°^ "^o ^'^'^ ^^^ refer to a pure power-law 
hypothesis (Tio)- The fit has been performed in the interval from 330 GeV-13 TeV. The last columns give the 
peak-emission energies -Emax- 



Set 


ON 

(h) 


4 


aPi 


(/'o 


a 


/3 


A 


-^Tnax 

(GeV) 


LF 

MF 
HF 


13.6 
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14.5 ±0.50 
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5 Source Intensity vs. Spectral Hardness 

Ne>e 
The hardness ratio is defined as i? = j^ siid. ^ where Ne^Ei^^ (resp. NE>En-,id) i^ the number of events 

^ low 

with a fitted energy greater than E'low (resp. E'mid)- 

All data have been divided into five sets with different average fluxes. For each data-set the hardness ratio 
has been computed for three different energy bands, [Eio^, E^idheV = {[450, 900]; [600, 1200]; [900, 1500]}, 
hereafter referred to as ^[>9oo/>450]' ^[>1200/>6001' ^nd ^[>i5oo/>900l- This method is more robust than the 
usual spectrum analysis for the following reasons: 

- since data have been divided into two bins in energy, statistical fluctuations on the single parameter 7? are 
reduced; 

- E'low is always chosen well-above the detector threshold in order to have a good estimation of the energy; 

- the data set is limited to zenith angles < 25° in order to avoid large threshold variations; 

- finally, an additional cut on the shower impact parameter (required to be < 130m) allows the energy resolu- 
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Figure 3: Hardness-ratio vs. source intensity for tliree different energy bands. Intensities are given as tlie integral flux above 
250 GeV in units of 10"" cm"^s"\ 



tion to be improved. 

The hardness ratio for the three different energy bands is shown in Figure 3, A correlation between hard- 
ness and intensity is clearly observed for all energy-bands. The x^ probabilities that the distributions are flat 



for -R[>900/>450] ' -^[>1200/>600] ^^l^ -R[>i500/>900] 



are: 3,8 x 10""^, 9.0 x 10~^, and 7 x 10"^ respectively. 



6 Conclusions 

The emission from Mrk 501 in 1997 and 1998 has been monitored with the CAT telescope. The spectrum 
has been measured from 330 GeV up to 10 TeV in 1997, revealing several spectral properties. A curved 
spectral shape has been derived, confirming the curvature reported by the Whipple group (Samuelson et al., 
1998) and HEGRA collaboration (Aharonian et al., 1999). The emission is seen to extend above 10 TeV. The 
observed intensity-hardness correlation can be simply ascribed to a shift of the peak TeV emission energy. 
This is confirmed by the observed increase in E'max derived from the spectral analysis. 
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